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MARKER AT THE ESTROGEN RECEPTOR GENE FOR 
DETERMINATI ON OP OSTEOPOROSIS PREDISPOSITION 

BACKGROUND OF THE INVENTION 

5 (a) Field of the Invention 

The invention relates to a method of determin- 
ing genetic predisposition to low or high bone mineral 
density of a patient, wherein low bone density is 
indicative of a predisposition to osteoporosis and high 
10 bone density is indicative of resistance to osteo- 
porosis . 

(b) Description of Prior Art 

Osteoporosis, a reduction of bone mineral den- 
sity (BHD), is a multifactorial disease that leads to 

15 an increasing risk of fracture and is becoming a major 
public health problem, especially in post-menopausal 
women. Its major consequence, hip fracture, has major 
health consequences with serious social, medical and 
economical implications in increasingly aging popula- 

20 tions. Research on osteoporosis emphasizes at either 
finding new therapeutic approaches or at the charac- 
terization of the major determinants of bone mineral 
density (BMD) in the hope to find markers useful in 
identification of women at risk of osteoporosis and its 

25 complications. Since therapy of established osteoporo- 
sis remains far from satisfactory, prevention is the 
best choice. 

Heredity has always been considered an impor- 
tant risk factor for osteoporosis, but its role was 

30 still poorly understood, at least until recently. 
Indeed, several studies of monozygotic (MZ) and dizy- 
gotic (D2) twins have shown that BMD was better corre- 
lated between MZ than DZ twins, indicating that BMD is 
genetically determined and that heredity could account 

35 for up to 80 to 90% of the variability in BMD. Intra 
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and intergeneration correlations were more contradic- 
tory, two studies showing significantly lower BMD in 
daughters of women with osteoporosis while another 
mother /daughter pairs study found no such difference. 
5 Since a decrease of one standard deviation in BMD 
within normal range approximately doubles the risk of 
fracture at different skeletal sites, the search for a 
marker likely to identify women at risk of post-meno- 
pausal osteoporotic fractures is becoming more and more 

10 relevant. 

More recently, a group of investigators have 
highlighted that polymorphism in the vitamin D receptor 
( VDR) gene were significantly correlated with the 
levels of osteocalcin, a marker of bone turnover. The 

15 same Australian group studied a cohort of 250 healthy 
twins for BMD and a polymorphism at the VDR revealed by 
the restriction enzyme Bsm I and found that genotype at 
the VDR could explain up to 75% of the genetic vari- 
ation in BMD at the lumbar vertebrae and proximal 

20 femur. The BB genotype, where B represents the allele 
with absence of the restriction site and b represents 
the allele with the polymorphic Bsm I site, is associ- 
ated with lower BMD and likely to an increased predis- 
position to osteoporosis. The bb genotype was associ- 

25 ated to a higher BMD and the Bb genotype was associated 
with intermediate BMD- Moreover, they studied 311 
unrelated healthy women and confirmed that genotype at 
the VDR locus was a strong predictor of BMD at both 
lumbar and femoral sites. The authors used an in vitro 
30 model (minigene model) suggesting that the VDR gene 
allelic variation could influence the rate of the 
receptor protein synthesis (Morrisson NA et al., 
Nature, 1994, 367:284-297). 

Many other groups have now published their 
35 results on differ nt populations which are contradic- 
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tory. A twin study done in the United Kingdom con- 
firmed the association between VDR genotype and BMD 
adjusted for body mass index (BMI). Another twin study 
in Indiana showed that up to 70% of the BMD could be 
5 attributed to heredity but did not find any relation- 
ship between VDR polymorphism and BMD. Another study 
concentrated on mother/daughter pairs in a Swiss popu- 
lation and found that VDR alleles contributed to 
mother/daughter BMD relationship. This study also 

10 showed that VDR polymorphism was associated with femo- 
ral BMD in teenage girls and pre-menopausal women. 
Another group from Massachusetts showed that the BB 
genotype was associated with significantly lower BMD 
compared to Bb and bb genotypes in both black and white 

15 races in a cohort of healthy unrelated pre-menopausal 
women aged 20 to 40. This suggested that the VDR poly- 
morphism may limit peak bone mass. 

Two studies from Japan on healthy unrelated 
women concluded that there was a difference in allele 

20 frequencies between Caucasian and Asian populations. 
One group found an association between VDR genotype and 
BMD at the lumbar vertebrae but the trend was more 
apparent in pre than post-menopausal women. The other 
* Japanese group concentrated on post-menopausal women 

25 and confirmed the association between VDR polymorphism 
and lumbar BMD but the genotype could not predict total 
BMD and had little clinical significance in the evalu- 
ation of bone status in Japanese women. Recently, a 
group from the Netherlands found an association between 

30 VDR alleles and BMD, but the allele previously 
described as "protective" against osteoporosis was 
associated with lower BMD at the femoral neck and 
Ward's triangle, suggesting allelic heterogeneity at 
the VDR locus. 
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Two studies looked at allelic frequencies in 
both normal and osteoporotic women. One group from 
Sweden found a 2.2 fold increase in the prevalence of 
the predisposing BB genotype in severely osteoporotic 
5 women compared to normal healthy controls but the rela- 
tionship did not seem to reach statistical signifi- 
cance. The other group from Minnesota, did not find a 
significant difference in genotype prevalence. They 
also noted that the effect of VDR genotype on BMD was 
10 modulated by age, with a greater effect in pre-meno- 
pausal women, reinforcing the hypothesis that allelic 
variations at the VDR locus may influence peak bone 
mass. 

Finally, a longitudinal study was recently pub- 
is lished where lumbar and proximal femoral BMD were meas- 
ured every six months for eighteen months in a cohort 
of elderly Swiss men and women. They concluded that 
the rate of bone loss was associated to VDR genotype at 
the lumbar spine but not at the femoral neck. BB 
20 genotype showed greater rate of bone loss p irrespective 
of calcium intake or supplement, but the rate of bone 
change for Bb genotype was influenced by overall 
calcium intake. 

The international Patent application No. 
25 WO 94/03633 (published on February 17, 1994) discloses 
a genetic test for assaying a predisposition to and/or 
resistance to high rates of bone turnover, development 
of low bone mass and responsiveness to therapeutic 
treatment. This test can be used for predicting osteo- 
30 porosis and likely response to preventive or therapeu- 
tic modalities- The test essentially consists in 
assessing the allelic variations in the vitamin D 
receptor gene. 
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SUMMARY OF THE INVENTION 

One aim of the present invention is to provide 
a method of determining predisposition to low or high 
bone density of a patient, which comprises determining 
5 estrogen receptor polymorphism in a biological sample 
of the patient, wherein heterozygosity is associated 
with high bone density and homozygosity is associated 
with low bone density. 

Another aim of the present invention is to pro- 
10 vide means to screen women to identify those for which 
the more expensive formal measurement of BMD is indi- 
cated. 

Another aim of the present invention is to pro- 
vide means of identifying young women that will be at 
15 risk of osteoporosis after their menopause so that they 
can attempt to increase their BMD to reach a higher 
peak bone mass. 

Another aim of the present invention is to pro- 
vide means of identification of target sub-groups of 
20 women for osteoporosis prevention measures/programs. 

Another aim of the present invention is to pro- 
vide means to determine which sub-group of post meno- 
pausal women will most benefit from osteoporosis treat- 
ments) and eventually predict their response to ther- 
25 apy or choose the optimal preventive pharmacotherapy. 

Another aim of the present invention is to 
identify means of prediction and management of BMD as 
well as biological parameters for the establishment of 
population-based osteoporosis prevention and interven- 
30 tion programs. 

In accordance with the present invention there 
is provided a method of determining predisposition to 
low or high bone mineral density and to development of 
osteoporosis of a patient, which comprises determining 
35 estrogen receptor polymorphism in a biological sample 
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of the pati nt, wherein heterozygosity is associated 
with high bone density and homozygosity is associated 
with low bone density. 

In some embodiments, the method of the present 
5 invention includes detecting the estrogen receptor 
polymorphism by analyzing the restriction fragment 
length polymorphism using an endonuclease digestion. 
The method can further include a step prior to the 
estrogen receptor gene digestion, wherein at least a 

10 fragment of the estrogen receptor is amplified, for 
example, by polymerase chain reaction. 

In accordance with the present invention, the 
estrogen receptor polymorphism, without limitation, is 
selected from the group consisting of PvuII polymorphic 

15 site located in the first intron of the ER gene or any 
DNA variant or mutation which shows some degree of 
linkage disequilibrium with one of the alleles of the 
PvuII polymorphism. 

The polymorphism of the estrogen receptor gene 

20 can be detected using at least one oligonucleotide spe- 
cific to the normal or variant estrogen receptor gene 
allele. 

In accordance with the present invention, low 
bone density is indicative of a predisposition to 
25 osteoporosis and/or bone fracture of the patient post- 
menopause, while high bone density is indicative of a 
resistance to osteoporosis and/or bone fracture of the 
patient post-menopause. 

In accordance with the present invention , 
30 estrogen receptor genotyping is indicative of response 
to therapy and/or to preventive treatments against low 
bone mineral density and bone and vertebrae fractures. 

The present invention also provides a kit for 
determining predisposition to low or high bone mineral 
35 density of a patient, which includes at least: 
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a probe sp cific for the estrogen receptor; 
an endonuclease selected from the group con- 
sisting of PvuII, PssI, Sad, and Xbal. 

Also, in accordance with the present invention, 
5 there is provided a method of determining predisposi- 
tion to estrogen hormone-related medical conditions of 
a patient, including the steps of: 

a) isolating nucleic acid of the patient from a 
biological sample; 
10 b) determining the genotype in said isolated 

nucleic acid of step a), wherein said genotype 
corresponds to a region including the gene 
encoding the estrogen receptor, and wherein 
heterozygosity or homozygosity is indicative of 
15 predisposition of said estrogen hormone-related 

medical conditions and homozygosity or het- 
erozygosity is indicative of likelihood of pro- 
tection against said estrogen hormone-related 
medical conditions. 
20 For the purpose of the present invention the 

following abbreviations and terms are defined below. 

The abbreviation "RFLP" refers to restriction 
fragment length polymorphism. 

The term "DNA polymorphism" or "polymorphic DNA 
25 sequence" refers to any sequence in the human genome 
that exists in more than one variant ( or version ) in 
the population. 

The term "estrogen hormone-related medical 
conditions" refers to, without limitation, any estro- 
30 gen-dependent diseases such as osteoporosis, endometri- 
osis, arteriosclerosis, breast cancer, ovarian cancer 
and other estrogen-dependent cancers, among others. 

The term "linkage disequilibrium" refers to any 
degree of non-random genetic association between one or 
35 more allele(s) of two different polymorphic DNA 
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sequences and that is du to the physical proximity of 
the two loci. Linkage disequilibrium is present when 
two DNA segments that are very close to each other on a 
given chromosome will tend to remain unseparated for 
5 several generations with the consequence that alleles 
of a DNA polymorphism (or marker) in one segment will 
show a non-random association with the alleles of a 
different DNA polymorphism (or marker) located in the 
other DNA segment nearby. Hence, testing of one of the 
10 DNA polymorphism (or marker) will give almost the same 
information as testing for the other one that is in 
linkage disequilibrium- This situation is encountered 
throughout all the human genome when two DNA polymor- 
phisms that are very close to each other are studied - 
15 Such a linkage disequilibrium with several 
polymorphisms in the vitamin D receptor gene are 
reported in Morrisson et al., 1994. Various degrees of 
linkage disequilibrium can be encountered between two 
genetic markers so that some are more closely associat- 
20 ed than others. 

The terms "estrogen receptor polymorphism" or 
"genetic marker" are intended to include, without limi- 
tation, PvuII (GDB (Genome Data Base) #G00-155-446) , 
Pssl (GDB #G00-155-447), SAcI (GDB #G00-155-448 ) , Xbal 
25 (GDB #G00-155-449) as well as the following ESR non- 
RFLP polymorphisms: GDB #G0O-162-450, #G00-162-541, and 
any other allelic variant of the estrogen receptor gene 
that shows some degree of linkage disequilibrium in any 
population sub-group with at least one of the above- 
30 mentioned estrogen receptor gene polymorphisms. 

The estrogen receptor gene polymorphism site in 
accordance with the present invention can be located 
within the estrogen receptor gene or within 500kb on 
each side of the ER gene, preferably within 250kb, and 
35 more preferably within 50kb. 
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BRIEF DESCRIPTION OF THB DRAWINGS 

Fig. 1A is a box plot graph of the ER/PvuII 

RFLP genotyping results correlated with bone mineral 
5 density at the L2-L4 vertebraes; and 

Fig- IB is a box plot graph of the ER/PvuII 

RFLP genotyping results correlated with bone mineral 
density at the hip. 

10 re-MMP DESCRIPTION OF THB IHVENTIOH 

The balance between bone formation and resorp- 
tion is very complex and is a consequence of many fac- 
tors (calcium and phosphor intake, PTH, estrogen's and 
androgens and 1, 25-dihydroxy vitamin D amongst others). 
15 Accelerated bone resorption is seen in the first few 
years after menopause characterized by ovarian failure 
and estrogen deficiency. We hypothesized that allelic 
variation at the estrogen receptor ( ER ) gene could be 
related to BMD. A cohort of 88 healthy unrelated post- 
20 , menopausal women aged between 60 and 70 were genotyped 
for a known PvuII DNA polymorphism in the first intron 
of the estrogen receptor (ER) gene (Yaich, Cancer Res., 
1992 r 52:77-83). A strong relationship between the ER 
genotype status and BMD was found, the heterozygotes 
25 for this RFLP having a significantly higher BMD than 
the two homozygote groups. This suggests that 
heterodimer rather than homodimer ER formation of this 
trans-acting element would confer a selective advantage 
against osteoporosis either through its interaction 
30 with the ligand and/or with DNA where it acts as a 
transcription factor. Furthermore, genotyping of the 
same subjects for a dinucleotide repeat (TA) n 5' to the 
ER locus (del Senno L et al., Hum. Mol. Genet., 1992, 
1(J>):354) showed a strong bimodal allelic distribution 
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and confirmed a decreased BMD in a subgroup of post- 
menopausal women but with less statistical power. 

These findings may have important repercussions 
on the screening for low BMD in post-menopausal women 
5 but also for the implementation of screening programs 
aimed at preventing low BMD in women and for identifi- 
cation of women which would benefit from estrogen 
replacement therapy at their menopause. Also, it may 
be useful for other applications in the prevention and 
10 treatment of osteoporosis but perhaps on other impor- 
tant diseases or conditions showing an estrogen-depend- 
ent behavior (cancer, arteriosclerosis and endometrio- 
sis amongst others). 

15 Subjects 

Eighty-eight (88) healthy unrelated post-meno- 
pausal women aged between 60 and 70 were recruited. 
They were all Caucasian French-Canadian women from a 
French background and were all living in the Quebec 

20 city metropolitan area. Recruitment was achieved 
through voluntary response to local newspaper adver- 
tisements for a study on factors affecting BMD in 
healthy post-menopausal women, including genetic fac- 
tors. The subjects had to answer a detailed question- 

25 naire covering family, medical, surgical, genealogical 
and obstetrical history. The survey also included 
information about medications and life habits 
(exercise, tobacco, alcohol, etc.). All women between 
60 to 70 who had accepted to sign consent form, answer 

30 adequately the questionnaire, have bone-mass measure- 
ment and a 30mL blood puncture were included unless 
they had a medical condition affecting bone homeosta- 
sis, or had used medications that modify bone metabo- 
lism. 
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Bon density 

BMD was measured at both the lumbar spine (L2- 
L4) and the femoral neck by dual-energy X-ray absorp- 
timetry (DEXA, DPX-L Lunar Radiation Corporation, Madi- 
5 son, Wisconsin, USA, Software version 3.2). This tech- 
nique is rapid, reliable and precise with a coefficient 
of variation of 1-2%. 

DMA isolation 

10 Blood samples were drawn into VACUTAINER™ con- 

taining EDTA and 200L was aliquoted in 1.5mL EPPENDORF 
™ tubes within 48 hours and stored at -20 °C until DNA 
purification. Genomic DNA was isolated from peripheral 
blood leukocytes with a minimethod necessitating only 

15 200L of whole blood where all steps are processed in a 
single 1.5mL tube as described below. Isolated DNA 
(5-7g) was resuspended into 100L TE 20:5 buffer (20mM 
Tris, 5mM EDTA), heated at 65°C for 4 hours and stored 
at 4°C. 

20 The whole procedure of DNA extraction takes 

place in a single 1.5mL EPPENDORF™ tube, minimizing 
sample identification errors and sample mixing. All 
reagent concentrations and volumes added to the tube 
are set in order to use a pipette repeat dispenser all 

25 through the procedure. The number of samples processed 
in a single DNA extraction/digestion procedure depends 
only on the availability of places in centrifuges 
turning at 13000 RPM. BIOFUGE™ 15 (Heraus) benchtop 
centrifuges in a 4°C room has a capacity of 80 x 1.5mL 

30 EPPENDORF™ tubes. The tubes remain in the 

centrifuge's 10-tube holders all through the procedure 
which eliminates a lot of tube manipulation. 

200nL of whole blood (collected on EDTA) is 
placed in an appropriately identified 1.5mL EPPENDORF™ 

35 tube. Th same tube will be used until gel loading so 
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it must be identified with a very resistant marking 
tool. A slightly heated metal tool may be used to 
engrave the number on the tube, this allows for perma- 
nent labeling. 

5 1 mL of a 20mM TRIS buffer + 5mM EDTA + 0,5% 

NP-40 solution <TE20:5+NP40) is added to the samples 
and they are left on ice for 30 minutes to allow for 
membrane disruption. Then the tubes are centrifuged at 
4000 RPM for 15 minutes and the supernatant is removed 
10 by gently shaking the tubes upside down over a sink. 
The pellets are then broken down by vortexing vigor- 
ously in a multitube vortex for a few minutes followed 
by two more washes with lmL TE20:5+NP40 and 15min. at 
4000 RPM followed by vortexing of the pellets at each 
15 time. Afterwards, the pellets are resuspended in lOO^iL 
TE20:5 (20mM TRIS buffer + 5mM EDTA) and 10^L of 10% 
sarcosine are added followed by 10\iL of 2mg/mL Pro- 
teinase K (BMC). The samples are incubated at 37°C 
overnight or 3 hours at 65°C. 
20 After the proteinase K digestion is completed 

100iiL ammonium acetate 7.5 M are added and well mixed. 
Then 500iiL of pure ethanol (cooled at -20°C) are added 
to each tube and mix well, allowing for DNA precipita- 
tion. There is enough DNA in each tube to actually see 
25 it precipitate when the cold ethanol is added. 

The tubes are centrifuged for 10 min. at 13 000 
RPM and the supernatant removed by gentle shaking over 
the sink. The pellet is then resuspended in 100 \iL 
TE20:5 (without NP40) and reprecipitated with 100 iiL 
30 ammonium acetate 7.5 M and 500 \iL of pure ethanol 
(cooled at -20°C). Tubes are centrifuged for 10 min. 
at 13 000 RPM and the supernatant again discarded- The 
DNA pellet is dried under negative pressure for 10 min. 
(in a desiccator), then resuspended in 200 \xh TE20:1 
35 (20mM TRIS buffer + ImM EDTA) and left at 37°C for a 
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few hours to allow for dissolution of the final DNA 
pellet. 

Pvu II polymorphism analysis 

5 From the resuspended DNA, 200ng of genomic DMA 

was amplified in a lOO^iL volume containing O.S^iM of 
both forward 5 • TGCCACCCTATCTGTATCTTTTCC3 1 (SEQ ID N0:1) 
and reverse 5 • TCTTTCTCTGCCACCCTGGCGTC (SEQ ID NO: 2) 
primers derived from Yaich ( Yaich, Cancer Res. , 1992, 

10 52:77-83), 200|iM of each of the four deoxyribonucleo- 
tides and 2.5U of Taq polymerase and buffer (Pr omega) 
in 1.5mM MgCl2- PCR amplification included the follow- 
ing steps: initial denaturation for 7 min. at 96°C 
followed by 35 cycles of amplification with denatura- 

15 tion at 94°C for 60s, annealing at 60°C for 60 s and 
polymerase extension at 72°C for 4 min. A final 
extension at 72°C for 10 min. was also included. 

About 500ng of the PCR product were digested 
with 10U Pvu II (New England Biolabs) overnight at 37°C 

20 and the polymorphism was visualized by ethidium bromide 
staining after a 2% agarose gel electrophoresis. 
Absence of the site (P) resulted into a 1.3Kb fragment 
whereas presence of the Pvu II site (p) resulted into 
850bp and 450bp fragments which was consistent with 

25 published data (Yaich, Cancer Res., 1992, 52:77-83). 
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Table 1 

DMA extraction protocol for miniprep of genomic DMA 

• Place 200nL of whole blood (collected on EDTA) in a 
1.5mL EPPENDORF 111 tube; 

5 • Add ImL of TE20:5 + NP40 ( 20mM TRIS buffer + 5mM EDTA 
+ 0,5% NP-40), leave 30 min. on ice and centrifuge 
for 15 min. at 4 000 RPM; 
•Throw away the supernatant; 

• Perform one more wash with ImL TE20:5+NP40, vortex 
10 and centrifuge for 15 min. at 4000 RPM; 

• Resuspend the pellet in 100 *iL TE20:5 and vortex 
well; 

• Add 10|iL of 10% sarcosine and then 10 \iL of Pro- 
teinase K(BMC) at 2.5 mg/mL and incubate at 37°C 

15 overnight or 3 hours at 65°C; 

• Add 100^L Ammonium acetate 7.5 M, mix well and then 
add 500^iL of pure ethanol (cooled at -20°C), mix well 
and centrifuge for 10 min. at 13 000 RPM; 

• Remove the supernatant and resuspend the pellet in 
20 100nL TE20:5 (without NP40); 

• Reprecipitate with 100 \lL Ammonium acetate 7.5M, mix 
well and 500jiL of pure ethanol (cooled at -20°C), mix 
well and centrifuge 10 min. at 13 000 RPM; 

• Remove the supernatant and dry the pellet by revers- 
25 ing the tubes (dry the cap of the EPPENDORF™ ) ; and 

• Detach the pellet in 400 \ih TE20:l(20mM TRIS buffer + 
ImM EDTA) . If the restriction enzyme to be used does 
not have a high temperature of digestion (50°C or 
more), incubate the tube at 65°C for a few hours 

30 before adding the restriction enzyme. 

Note: It is possible to use the same 1.5mL EPPENDORF™ 
tube from the beginning to the end of the procedure. 
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Mlcr satellite <TA) n polymorphism analysis 

The 88 subjects were subsequently genotyped for 
a known dinucleotide microsatellite (TA) n 5 1 to the ER 
gene (GDB # G00-162-541) that had 17 alleles and showed 
5 82% heterozygosity in a previous report. PGR was car- 
ried out in an MJ PTC-100™ thermocycler with hot-bon- 
net (MJ Research Inc, Watertown, Ma). Each 50^iL reac- 
tion contained 250 ng of genomic DNA, 40nM of end- 

32 

labeled with [gamma PldATP forward primer 

10 5' GACGCATGATATACTTCACC 3' (SEQ ID NOr3) and 200nM of 
reverse primer 5* GCAGAATCAAATATCCAGATG 3' (SEQ ID 
N0:4), 200tiM of each dNTP, 3.5mM MgCl2, 15% glycerol 
and 2.5U of ULTRATHERM™ DNA polymerase and buffer 
(Biocan Scientific Inc, Mississauga, Ont). Amplifica- 

15 tion conditions were: initial denaturation for 7 min. 
at 96 °C followed by 30 cycles of amplification with 
denaturation at 94°C for 30s, annealing at 50°C for 45s 
and polymerase extension at 72 °C for 60s. A final 
extension at 72°C for 2 min. was included. 5^L of the 

20 PCR reaction were then mixed in the same volume of 
loading buffer (Bromophenol blue) and 5-6|iL of the 
mixture were deposited and the radiolabeled products 
were resolved on a 6% denaturing polyacrylamide gel 
electrophoresis, exposed overnight at room temperature 

25 and visualized by autoradiography. 

Statistical analysis 

ANOVA analyses were performed given the fact 
that the dependent variables were continuous (bone den- 

30 sity measurements at different sites) and the inde- 
pendent variables were ordinal (genotypes at the VDR 
and ER genes). The significance level was set at alpha= 
0.05. A Chi square analysis was performed to study the 
level of linkage disequilibrium between the two differ- 

35 ent ER polymorphisms. 
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Results 
PvuII RFLP 

Typing of the ER receptor gene alleles using PCR 
5 followed by PvuII digestion of the PCR products was 
performed on the 88 samples of women between 60 and 70 
years old. This RFLP has two alleles, namely p and P 
and each individual having two chromosomes 6 carries 
two copies of the ER gene and hence , two alleles of the 

10 RFLP. The various combinations of these two possible 
alleles generates three different genotypes at this 
locus: pp, pP and PP. Where "pp" designates homozy- 
gotes for the presence of the PvuII site, H PP" is for 
homozygotes for the absence of the PvuII site, and "pP" 

15 is for heterozygotes for the PvuII site with one allele 
with the presence of the site and the other with the 
absence of the site. 

ER/PvuII RFLP genotyping results were correlated 
with bone mineral density as determined by DEXA at the 

20 L2-L4 vertebrae (Fig. 1A) as well as at the hip (Fig. 
IB). ANOVA analysis showed that the ER genotype was a 
strong predictor of BMD at the L2-L4 level (p=0.0097) 
composed mainly of trabecular bone as well as at the 
hip (p=0.0387), which is composed principally of 

25 cortical bone. In fact, the ER genotype could divide 
women in two groups with a mean difference in BMD at 
L2-L4 of 0.12 g/cm3 (11%) and at the hip of 0.08 g/cm3 
(10%). This represents a one standard deviation 
difference between the mean BMD at each site between 

30 the two groups of women as specified by the ER genotyp- 
ing. 

Unexpectedly , individuals homozygous for either 
of the ER alleles (pp or PP) had comparable mean BMD 
(no statistical difference, p>0.05) at both sites stud- 
35 ied. However, their BMD was one standard deviation 
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below those women heterozygous for the ER g nes (pP). A 
one standard deviation difference in bone mineral den- 
sity in elderly women is well known to be associated 
with an increased relative risk of bone fracture of 
5 about 2.4. 

Hence, ER genotyping using the PvuII RFLP system 
allows to classify women in two groups (homozygotes vs. 
heterozygotes ) that have a more than two-fold dif- 
ference in the risk of bone fracture - 

10 

Microsatellite 

In order to test for the hypothesis of a founder 
effect in the variation at the ER genotype observed in 
the population, all samples were retested with a ER 

15 gene microsatellite. This highly polymorphic genetic 
marker had a bimodal allele distribution (Table 2) and 
hence the alleles were grouped into two categories 
(E=alleles D to G? M = alleles H to P) to form a bial- 
lelic system generating three possible genotypes (EE, 

20 EM, MM) . The microsatellite genotypes were highly cor- 
related with the RFLP genotypes (Table 3) (p<0.0001) 
confirming a limited number of ER alleles present in 
the population. However, the linkage disequilibrium 
observed was not absolute (Table 2). When the micro- 

25 satellite genotypes were correlated with BMD measure- 
ments, only a trend was observed (p=0.11), one of the 
between group comparisons was barely significant (EM 
vs. MM: p=0. 05, mean difference in BMD=0.09 g/cm2) but 
in the same direction as the strong correlation 

30 observed with the PvuII RFLP. 
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Table 2 

Allelic frequencies for BR -MS 
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Absolute 


10 


45 


10 
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13 
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33 
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12 
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(%) 





























5 Trtle 3 

Contingency table of PvuII RPLP vs ER -MS 



Observed Frequencies for ER, ERms 



EE EM MM Totals 

pp 20 6 2 28 

pp 3 25 4 32 

pp 2 2 19 23 

Totals 25 33 25 83 



10 There is shown in Table 3 an example of the 

linkage disequilibrium between the alleles p and E, 
where the homozygotes EE are associated with the 
homozygotes pp. 

15 Discussion 

Previous studies have reported the effect of 
genetic variation of the VDR gene on BMD in women and 
the results appear contradictory. Whether these dis- 
crepancies are due to differences in the samples stud- 

20 ied, or differences in methodological or data analysis 
approaches, or correspond to real differences between 
the populations studied has to be resolved. We did not 
find any effect of the VDR genotype on BMD of those 
women included in the present study but the effects of 

25 genetic variation of another important steroid receptor 
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gene, namely the estrogen receptor (ER) , were also ana- 
lyzed . 

In accordance with the present invention, it is 
demonstrated that genotyping at the ER gene allows one 
5 to predict whether women will be at increased risk of 
osteoporosis and bone fracture when they reach 60 years 
old. Because the ER genotype is genetically determined 
and remains the same throughout life, it is possible to 
genotype at the ER gene young women who have not yet 

10 reached their peak bone mass and concentrate preventive 
actions to increase BMD or diminish bone loss for those 
women which have a homozygous ER genotype as they will 
have, as a group, a BMD one S.D. below heterozygotes 
(i.e. more than twice the risk of osteoporosis) when 

15 they reach 60 years of age. 

It is also important that the ER genotype was 
correlated with both BMD of trabecular bone and corti- 
cal bone which are two metabolically different types of 
bone. Furthermore, the intensity of the effect of the 

20 ER genotype was similar (i.e. one Standard Deviation) 
for each type of bone. ER genotyping thus may be a 
good marker of the homeostatic set point of general 
bone metabolism. 

It is worth mentioning the unexpected relation- 

25 ship between the ER genotypes and BMD at the hip and 
vertebrae. The results reported by Morrisson et al. on 
the relationship of VDR genotypes with BMD showed a 
cumulative (or dosage) effect of one allele over the 
other i.e. that the BMD increased with the number of 

30 copies of a given allelic variant of the VDR gene 
(Morrisson NA et al., Nature, 1994, 367:284-297). The 
present analyses of the ER genotypes revealed an inter- 
action (or inverted- M v" shaped) effect instead of a 
cumulative effect with the best predictor of a low BMD 

35 being the presence of identical (homozygous) alleles of 
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the ER gene. This group of women (pp or PP) repre- 
sented 60% of the sample studied. 

This simple genetic test of the ER genotype 
could also potentially be used to identify, prior to 
5 their menopause or later, which women would most bene- 
fit of estrogen replacement therapy or preventive phar- 
macotherapy ( such as with biphosphonates ) . Further 
studies may also reveal that the clinical response to 
such therapies could be related to the ER genotype. 

10 Also, ER geno typing, as it is inexpensive, may be used 
to screen post-menopausal women to identify a sub-group 
which is at higher risk of low BMD and who could 
benefit from a more formal BMD measurement such as DEXA 
(which is too expensive to be offered as a screening 

15 procedure). Hence genotyping of the ER could become a 
fundamental parameter in prediction and management of 
low BMD as well as for the establishment of population- 
based osteoporosis prevention and intervention 
programs . 

20 Also, typing of both the VDR and ER gene poly- 

morphisms may enable a better prediction of BMD even if 
we did not find such a performance in the population 
studied. It is however likely that the association of 
ER genotyping with other analytical procedures 

25 {measurement of bone metabolites, other genotypes, 
etc.) may allow an even better discrimination between 
women of high and low risk for osteoporosis - 

The biological mechanisms by which this strong 
effect of the ER genotype on trabecular and cortical 

30 BMD take place remains unknown. However, given the 
very peculiar correlation between the genotypes and the 
BMD measurements, it is tempting to speculate that the 
lower BMD in ER homozygotes vs ER heterozygotes may 
indicate that there is a difference in the physiologi- 

35 cal performance of heteroduplex estrogen receptors as 
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compared to homoduplex estrogen receptors. One possi- 
ble mechanism is, hence, that the mutation (or polymor- 
phism) involved affects the dimerization domains of 
receptor monomers and influences the control of bone 
5 metabolism by the steroid hormone. This could be con- 
firmed once the mutation involved in this ER effect on 
BMD after 60 years of age is identified; probably, the 
PvuII RFLP is closely associated with the effective ER 
gene mutation that has yet to be identified by sequenc- 
10 ing. 

Other mechanisms of action of ER gene polymor- 
phism can also be postulated, including the effect of 
heterodimeric receptors on hormone or DNA binding or 
even on binding with other proteins involved in the 

15 availability or efficiency of ER receptors for the hor- 
monal control of genes and cellular processes. 

The present work demonstrating ER genotype 
effects on a disease clearly associated with estrogen 
metabolism opens the field of other estrogen-dependent 

20 diseases/conditions such as arteriosclerosis , 
endometriosis, breast cancer, ovarian cancer, and other 
estrogen-dependent cancers. 

We have demonstrated for the first time that 
estrogen-dependent cellular processes can vary from one 

25 individual to the other according to the combination of 
estrogen receptor variants (or to. the estrogen receptor 
genotype). Since there is only one estrogen receptor 
gene per haploid human genome (each cell contains two 
haploid genomes ) , the differences in biological 

30 efficacy of the estrogen hormone via the estrogen 
receptor genotype disclosed in the present application 
will also apply to other estrogen dependent diseases or 
conditions. 

While the invention has been described in con- 
35 nection with specific embodiments thereof, it will be 
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understood that it is capable of further modifications 
and this application is intended to cover any varia- 
tions, uses, or adaptations of the invention following, 
in general, the principles of the invention and 
5 including such departures from the present disclosure 
as come within known or customary practice within the 
art to which the invention pertains and as may be 
applied to the essential features hereinbefore set 
forth, and as follows in the scope of the appended 
10 claims. 
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SEQUENCE LISTING 



(1) GENERAL INFORMATION: 

(i) APPLICANT: 

(A) NAME: UNIVERSITE LAVAL 

(B) STREET: Cite Universitaire 

(C) CITY: Quebec 

(D) STATE: Quebec 

(E) COUNTRY: Canada 

(F) POSTAL CODE (ZIP) : G1K 7P4 

(G) TELEPHONE: (418) 656-3525 

(H) TELEFAX: (418) 656-7785 

(A) NAME: ROUSSEAU f Francois 

(B) STREET: 2540 avenue Lalande 

(C) CITY: Ste-Foy 

(D) STATE: Quebec 

(E) COUNTRY: Canada 

(F) POSTAL CODE (ZIP): G1W 1M7 



(ii) TITLE OF INVENTION: MARKER AT THE ESTROGEN RECEPTOR GENE 
FOR DETERMINATION OF OSTEOPOROSIS PREDISPOSITION 

(ill) NUMBER OF SEQUENCES: 4 

(v) COMPUTER READABLE FORM: 

(A) MEDIUM TYPE: Floppy disk 

(B) COMPUTER: IBM PC compatible 

(C) OPERATING SYSTEM: PC-DOS /MS -DOS 

(D) SOFTWARE: Patentln Release #1.0, Version #1.30 (EPO) 

(vi) PRIOR APPLICATION DATA: 

(A) APPLICATION NUMBER: US 08/592,835 

(B) FILING DATE: 25-JAN-1996 



(2) INFORMATION FOR SEQ ID NO:l: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 24 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO:l: 
TGCCACCCTA TCTGTATCTT TTCC 24 



(2) INFORMATION FOR SEQ ID NO:2: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 23 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 
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(D) TOPOLOGY; linear 
(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 2: 
TCTTTCTCTG CCACCCTGGC GTC 



(2) INFORMATION FOR SEQ ID NO: 3: 

(i) SEQUENCE CHARACTERISTICS: 
{A) LENGTH: 20 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 3: 
GACGCATGAT ATACTTCACC 



(2) INFORMATION FOR SEQ ID NO: 4: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 21 base pairs 

(B) TYPE: nucleic acid 
<C) STRANDEDNESS: single 
(D) TOPOLOGY: linear 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 4: 
GCAGAATCAA ATATCCAGAT G 
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WHAT IS CLAIMED IS ; 

1. A method of determining predisposition to low 
or high bone density of a patient, which comprises 
determining estrogen receptor polymorphism in linkage 
disequilibrium in a biological sample of said patient, 
wherein heterozygosity is associated with high bone 
density and homozygosity is associated with low bone 
density. 

2. The method of claim 1, wherein detecting said 
estrogen receptor polymorphism comprises analysis of a 
restriction fragment length polymorphism using endonu- 
clease digestion. 

3. The method of claim 2, which further comprises 
a step prior to said estrogen receptor gene digestion, 
wherein at least a fragment of said estrogen receptor 
is amplified. 

4. The method of claim 3, wherein said estrogen 
receptor polymorphism is selected from the group con- 
sisting of PvuII polymorphic site located in the first 
intron of the ER gene or any DNA variant or mutation 
which shows some degree of linkage disequilibrium with 
one of the alleles of the PvuII polymorphism. 

5. The method of claim 4, wherein said polymor- 
phism of the estrogen receptor gene is detected using 
at least one oligonucleotide specific to the normal or 
variant estrogen receptor gene allele. 

6. The method of claim 3, wherein said fragment of 
said estrogen receptor is amplified by polymerase chain 
reaction. 
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7. The method of claim 1, wherein low bone density 

is predisposition to osteoporosis and/or bone fracture 
of said patient during post-menopause. 

8* The method of claim 1, wherein high bone den- 

sity is indicative of resistance to osteoporosis and/or 
bone fracture of said patient during post-menopause. 

9. The method of claim 1, wherein estrogen recep- 
tor genotyping is indicative of response to therapy 
and/or to preventive treatments against low bone min- 
eral density and bone and vertebrae fractures. 

10. A prognosis kit for determining predisposition 
to low or high bone mineral density of a patient, which 
comprises at least a probe specific for estrogen recep- 
tor; an endonuclease selected from the group consisting 
of PvuII, PssI, Sad, and Xbal. 
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